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SECTION  1.0 


INTRODUCTION 


1.  1  Objective 

The  objective  of  this  program  was  to  evaluate  the  capability  of  the  NET-2 
Circuit/System  Analysis  Computer  Program^  as  an  aid  in  the  analysis 
of  radiation  effects  on  complex  semiconductor  devices  and  microcircuits. 
Mathematical  models  considered  included  both  the  terminal  built-in  models 
and  lumped  models  of  bipolar  and  MOS  devices.  The  evaluation  included 
the  definition  of  each  model  in  terms  acceptable  to  the  CDC  6600  version 
of  NET-2,  computation  of  the  electrical  performance  and  transient  radiation- 
induced  response,  comparison  of  the  calculated  results  to  "exact"  sohitions 
when  possible  to  determine  accuracy,  and  comparison  of  program  execution 
time  to  calculations  with  other  models  to  judge  the  computational  efficiency. 

1 . 2  Technical  Approach 

There  are  many  semiconductor  device  analysis  techniques,  models,  and 
computer  programs  to  aid  in  the  numerical  analysis.  The  device  models 
and  elements  provided  in  NET-2  allow  the  convenient  use  of  greater  detail 
and  accuracy  in  simulation  than  that  allowed  in  other  circuit  analysis  pro¬ 
grams  (e.g. ,  CIRCUS,  SCEPTRE,  SYSCAP).  This  capability  is  not  as 
complete  as  that  available  in  special-purpose  computer  programs  (such  as 
those  developed  by  BTL  and  IBM),  but,  unlike  the  device  analysis  programs, 
enables  simulation  of  the  device  under  fully  general  circuit  conditions. 

In  the  evaluation  of  the  NET-2  semiconductor  device  analysis  capability  it 
was  first  necessary  to  develop  confidence  in  both  NET-2  and  the  lumped 
model  elements  necessary  for  the  generalized  device  modeling.  This  was 
accomplished  by  first  considering  the  electrical  performance  and  radiation- 
induced  photoresponse  of  simple  devices  using  familiar  first-order  mathe¬ 
matical  models.  The  lumped  model  representations  of  the  device  could  then 
be  compared  directly  to  familiar  results.  The  next  step  was  to  increase  the 
complexity  of  the  simple  device  models  by  using  more  lumped  model  elements 
and  including  effects  (such  as  built-in  electric  fields)  that  are  generally  be¬ 
yond  the  scope  of  the  familiar  terminal  device  model.  Once  the  credibility 
of  the  NET-2 /lumped  model  technique  had  been  established,  study  was  con¬ 
tinued  to  establish  the  perspective  necessary  to  suggest  rewarding  applica¬ 
tions  of  NET-2  to  existing  analytical  techniques.  This  included  the  derivation 
of  complex  models  for  the  elements  of  a  junction-isolated  bipolar  mieiocircuit 
(including  the  multiple-emitter  transistor)  and  the  analysis  of  a  complete 
junction-isolated  TTL  Gate  microcircuit. 

T  A.  F.  Malmberg,  "NET-2  Network  Analysis  Program  -  Preliminary 
User's  Manual,  ”  Harry  Diamond  Laboratories;  May  1970. 
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1.  3 


Summary  of  Results 


All  the  lumped  model  elements  (combinance,  storance,  diffusance,  driftance, 
and  carrier  generation)  can  be  used  with  the  built-in  p-n  junction  model  to 
simulate  the  electrical  performance,  transient  ionizing  radiation-induced 
photoresponse,  and  time-dependent  neutron-induced  displacement  damage 
effects  in  virtually  any  bipolar  semiconductor  device.  In  the  present  study, 
device  operation  under  conditions  of  junction  voltage  breakdown  were  not 
considered  and  some  modifications  would  be  necessary  in  the  p-n  junction 
model.  The  terminal  device  models  evaluated  were  those  for  the  p-n  junction 
diode,  the  zener  diode,  the  bipolar  transistor,  and  the  MOS  transistor,  all 
found  to  be  computable  and  accurate  to  the  limitations  implicit  within  the  model. 

Overall  impressions  on  the  use  of  NET-2  were  very  favorable.  Initial  con¬ 
vergence,  numerical  stability,  accuracy  of  the  solution  and  computer  running 
times  were  all  good.  Computer  time  for  a  "typical"  transient  response  cal¬ 
culation,  d-c  characterization,  or  swept  frequency  response  were  on  the 
order  of  100  seconds  for  problems  involving  substantial  complexity  to  the 
analyst,  but  modest  in  terms  of  the  code's  capacity  (e.  g.  ,  25  nodes,  100 
elements).  The  cost  trade-off  in  using  the  lumped  model  representation 
rather  than  a  terminal  model  appeared  to  be  somewhat  independent  of  the 
number  of  lumped  model  elements  used  and  more  dependent  on  the  number 
of  p-n  junction  models  used.  It  appeared  that  each  p-n  junction  model  and 
associated  lumped  model  elements  was  roughly  equivalent  to  a  single  ter¬ 
minal  transistor  model.  The  accuracy  of  the  lumped  model  can,  however, 
often  be  much  greater  than  a  factor  of  two  better  than  that  obtained  with  a 
terminal  transistor  model. 


SECTION  2.  0 

LUMPED-MODEL  TECHNIQUE 


The  lumped  model  technique  is  simply  the  application  of  a  difference-differ¬ 
ential  approximation  to  the  partial  differential  equations  describing  carrier 
distribution  and  motion  in  a  semiconductor.  Given  a  particular  semiconductor 
device,  each  bulk  semiconductor  region  is  characterized  in  terms  of  a  "net¬ 
work”  representing  carrier  distribution  and  current  throughout  the  region. 

The  p-n  junction  models  provide  boundary  conditions  for  each  region  relating 
the  earner  density  at  the  edges  of  the  junction  to  the  applied  junction  voltage. 
The  general  lumped  model  technique  proposed  by  Linvill^*  3.  4  includes  rep¬ 
resentation  of  both  minority  and  majority  carrier  flow  by  either  signal- 
dependent  diffusion  or  drift.  The  lumped  model  representation  of  NET-2 
has  been  simplified  to  gain  more  efficient  computability  and  allows  charac¬ 
terization  of  the  semiconductor  region  by  the  minority  carrier  distribution 
and  flow.  This  is  essentially  the  assumption  of  space-charge  neutrality  and 
is  implicit  in  virtually  all  bipolar  device  models.  ^ 

The  principal  advantages  of  the  NET-2  /lumped  model  are:  1)  direct  corre¬ 
spondence  to  the  physical  processes  and  radiation  effects  of  the  bulk  semi¬ 
conductor,  2)  flexibility  in  the  definition  of  the  model  complexity,  and  3) 
compatibility  in  the  engineering  analysis  of  many  interconnected  devices, 
circuits,  or  system  blocks. 

The  elements  and  parameters  of  the  lumped  model  are  illustrated  in  Figure 

1  for  a  one-dimensional  representation  of  a  bulk  semiconductor  region. 
Carrier  recombination,  storage  and  generation  in  each  sub-region  (lump) 
are  represented  by  the  lumped  model  combinance  (Hc),  storance  (Sj)  and 
current  generator  (i^)  elements,  respectively.  Current  flow  between  lumps 
by  carrier  diffusion  and/or  drift  is  represented  by  the  current  through  the 
diffusance  (H(j)  and  driftance  (Df)  elements  respectively.  It  is  important  to 
note  that  the  current  flow  in  all  elements  is  defined  as  positive  in  the  direction 
of  "conventional"  current  flow.  When  electrons  are  the  minority  carriers, 
the  flow  of  carriers  is  then  opposite  the  flow  of  conventional  current. 

It  is  helpful  to  first  just  consider  the  carrier  generation  and  recombination 
processes  in  n-  and  p-type  semiconductor  as  sh*>wn  schematically  in  Figure 
2.  The  minority  carrier  density  is  specified  in  terms  of  the  "excess"  density, 
or  that  value  different  from  the  thermal  equilibrium.  For  the  k-th  lumped 
region  then, 

2  J.  G.  Linvill,  "Lumped  Models  of  Transistors  and  Diodes,  "  Proc.  IRE, 
Vol.  46,  pp.  1141-1152;  June  1958. 

J.  G.  Linvill  and  J.  F.  Gibbons,  Transistors  and  Active  Circuits,  New 
York:  McGraw-Hill  Co. ,  1961. 

^  J.  G.  Linvill,  Models  of  Transistors  and  Diodes,  McGraw-Hill,  New 
York;  1963. 

^  J.  P.  Raymond  and  R.  E.  Johnson,  "Detailed  Lumped-Model  Analysis  of 
Transistor  Ionizing  Radiation  Effects,"  IEEE  Trans.  Nucl.  Sci.  ,  Vol. 

NS- 13,  No.  6  pp.  95-104;  December  1966. 
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k  th  lump  (sub -region) 


(a)  One-di  mensional  semiconductor  region 


Dfk 

(b)  NET-2  lumped  model  elements:  n-type  semiconductor 


(c)  NET-2  lumped  model  elements:  p-type  semiconductor 


Figure  1  One-Dimensional  Lumped  Model 
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p,  =  P,  -  p  for  n-tyoe  semiconductor 
k  k  n 

n  =  N  -  n  for  p-type  semiconductor 

K  K  p 

where  p^  and  n^  are  the  excess  minority  carrier  densities,  and  are  the 
actual  minority  carrier  densities,  and  pn  and  iip  are  the  thermal  equilibrium 
densities.  The  process  of  a  single  event  carrier  recombination  in  a  semicon¬ 
ductor  removes  one  hole  and  one  electron  as  free  carriers.  If  the  electron 
density  is  reconsidered  as  a  negative  charge,  then  the  process  of  recombina¬ 
tion  is  the  transfer  of  a  positive  charge  from  the  hole  side  of  the  combinance 
to  the  electron  side  of  the  combinance  as  shown  in  Figure  2.  This  reduces 
the  absolute  value  of  both  the  hole  density  (a  positive  number  subtracted 
from  a  positive  quantity)  and  the  electron  density  (a  positive  number  added 
to  a  negative  quantity).  The  carrier  recombination  rate  then  is  a  conventional 
current  flowing  from  the  hole  side  of  the  lumped  model  to  the  electron  side 
and  is  given  by, 


and 


i  =  p  H  for  n-type  semiconductor 
r  k  cp 

P 


H  =  qAMV/x 
cp  p 


1  =  n.  H  for  p-t.ype  semiconductor 

r  k  cn 

n 


H 

cn 


qAiW/T 

n 


Physically,  this  simply  represents  the  average  carrier  recombination  in  the 
lump  as  proportional  to  the  excess  density,  proportional  to  the  volume  (Ai-W), 
and  inversely  proportional  to  the  minority  carrier  lifetime  (t).  Neutron- 
induced  permanent  damage  effects  are  included  in  the  lumped  model  by  direct 
degradation  of  the  minority  carrier  lifetime  or, 

-r—  (t)  =  ~ -  +  f  [  f  (t),  K(t)  1 
f  i 


where  "f  is  the  resultant  lifetime,  r.  is  the  initial  lifetime,  !  (t)  is  the  time- 
dependent  neutron  fluence,  and  K{t)  is  the  lifetime  damage  constant  (time- 
dependent  due  to  short-term  annealing  effects). 

Ionizing  radiation  effects  are  represented  in  the  lumped  model  directly  by 
including  the  radiation-induced  carrier  generation  in  each  lump.  For  high- 
energy  radiation,  the  ionization  will  be  proportional  to  the  radiation  intensity, 
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and  will  be  absorbed  uniformly  throughout  the  device.  Since  the  process  of 
carrier  generation  is  opposite  that  of  carrier  recombination,  the  current 
generators  of  the  lumped  model  are  in  the  direction  opposite  that  of  the 
carrier  recombination  current,  as  shown  in  Figure  2.  Thus,  for  both  n- 
type  and  p-type  material  we  have, 

i  =  qAAWg  V  (t) 
gk  o 

where  y  (t)  is  the  time-dependent  absorbed  dose  rate. 

In  addition  to  representing  carrier  generation  and  recombination  in  each 
lumped  region,  we  must  account  for  the  increase  or  decrease  of  carriers 
under  transient  conditions.  This  is  represented  in  the  lumped  model  by 
storance  element  (Sj.  and  Stn  in  Figure  2),  which  acts  much  like  a  capaci¬ 
tance  for  minority  carriers.  The  current  through  the  storance  is  given  by, 

for  n-type  semiconductor 

for  p-type  semiconductor 

Carrier  transport  between  lumps  is  represented  by  the  carrier  diffusance, 
for  carrier  transport  by  diffusion,  and  carrier  driftance,  for  carrier  trans¬ 
port  as  the  result  of  a  built-in  electric  field.  The  current  through  the  dif¬ 
fusance  element  is  given  by, 


where 


H 


=  (Pl.  '  Pl- 


k  + 


1>  Hdk 


H  =  qAD/MV 


and  D  -  minority  carrier  diffusion  constant 

The  direction  of  current  flow  in  the  diffusion  is  shown  in  Figure  3  for  n- 
type  semiconductor  and  p-type  semiconductor.  When  the  minority  carriers 
are  electrons,  the  direction  of  conventional  current  flow  will  be  in  the  direc¬ 
tion  opposite  that  of  the  electron  diffusion. 


The  current  through  the  driftance  clement  is  given  by, 

1. 


D 


fk 


,pk  t  i  1  pk>  Dtk 


where 


D 


fk 


qA  u  E 
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u  =  minority  carrier  mobility 

E  =  built-in  electric  field,  positive  in  the  direction 
of  the  driftance  arrow 

It  is  important  to  note  that  the  driftance  element  is  oriented  in  the  direction 
of  positive  electric  field  in  both  n-type  and  p-type  semiconductor.  For  the 
p-type  semiconductor  then,  the  direction  of  current  through  the  driftance 
will  be  opposite  to  that  of  electron  flow,  and  the  orientation  of  the  driftance 
will  be  opposite  to  the  direction  of  positive  acceleration  for  the  electrons. 
Current  flows  through  the  driftance  elements  in  n-type  and  p-type  semicon¬ 
ductor  are  shown  in  Figure  3. 

All  the  lumped  model  elements  available  in  NET-2  are  summarized  in  Fig¬ 
ure  4.  The  number  of  lumps  used  to  characterize  a  bulk  semiconductor  region 
of  a  bipolar  device  varies  from  1  to  about  10,  depending  on  the  accuracy  de¬ 
sired  and  the  degree  of  variation  of  built-in  electric  field  and  bulk  semicon¬ 
ductor  parameters  through  the  region.  Representations  of  a  single  region 
by  a  model  containing  more  than  10  lumps  is  probably  not  realistic  because 
of  the  NET-2  numerical  accuracy  limitations. 


Models  of  each  bulk  semiconductor  region  of  the  device  are  connected  to  each 
other  with  the  built-in  p-n  junction  model,  or  are  connected  to  the  external 
circuit  terminals  with  an  infinite  combinance  (shor ..-circuit)  representing  an 
ohmic  contact.  The  p-n  junction  model  of  NET-2  is  shown  in  Figure  5.  The 
minority  carrier  excess  density  at  the  edges  of, the  junction  is  defined  in  terms 
of  the  applied  junction  voltage, 


p  .  exp  (•'  V )  -  1  ] 
n 


for  the  n-side 


n 


1 


n  "  exp  (vV)  -  1  ' 

p 


for  the  p-side 


where  v  =  q/kT  and  pn,  n  arc  the  thermal  equilibrium  minority  carrier  den¬ 
sities  at  the  n-type  and  p-type  sides  of  the  junction  respectively.  Currents 
through  the  junction  are  continuous.  That  is,  the  current  that  enters  one 
side  as  a  minority  carrier  flow  is  identical  to  that,  leaving  the  other  side  of 
the  junction  as  a  majority  carrier  flow.  The  depletion  capacitance  of  the 
junction  is  connected  between  majority  carrier  nodes.  Doping  levels  and 


grading  of  the  junction  determine  the  zero  volt  capacitance  (C0),  the  grade 
constant  n  (1/2  for  an  abrupt  junction  to  1/3  for  a  linearly  graded  junction) 


and  barrier  potential,  V  . 


The  lumped  model  elements  and  p-n  junction  model  described  enable  the 
NET-2  user  to  define  a  wide  range  of  models  for  many  bipolar  devices  and 
elements  of  monolithic  integrated  circuits.  The  following  discussion 
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Figure  4.  Summary  of  Lumped-Model  Parameters 
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Figure  5.  NET-2  P-N  Junction  Model 


(Section  3)  presents  a  summary  of  the  results  that  lead  to  the  credibility  of 
NET-2 /lumped  model  analysis.  Section  4  presents  the  detailed  format  (with 
comments)  for  two  semiconductor  device  analysis  problems  that  may  be 
helpful  to  the  user  as  examples. 


12 


SECTION  3.  0 


TECHNICAL  RESULTS" 

3.  1  P-N  Junction  Diode  Analysis 

J 

Electrical  behavior  and  photoresponsc  of  an  ideal  p'  -n  diode  can  be  charac¬ 
terized  by  determining  the  carrier  flow  and  distribution  in  the  high-resistivity 
n-region  of  the  diode.  The  general  form  of  the  lumped  model  diode  represen¬ 
tation  is  shown  in  Figure  t>.  Combinance  and  storancc  elements  represent 
the  average  of  carrier  recombination  and  storage  for  a  finite  region  of  bulk 
semiconductor.  Complex  ionizing  radiation  effects,  as  observed  from  the 
overall  device  photoresponsc,  are  represented  simply  as  the  radiation-induced 
carrier  generation  proportional  to  the  time-dependent  ionizing  radiation  ab¬ 
sorbed  dose  rate.  Carrier  transport  between  lumped  regions  by  either  dif¬ 
fusion  or  drift  is  represented  in  the  lumped  model  by  the  current  through  the 
diffusancc  and  driftance  elements,  respectively.  The  lumped  model  repre¬ 
sentation  of  the  diode  n-region  is  defined  on  one  side  by  the  p-n  junction 
model,  and  on  the  other  side  by  the  ohmic  contact  (or  infinite  recombination 
surface).  The  number  of  lumps  used  in  the  representation  and  their  relative 
geometry  is  arbitrary  and  influences  the  accuracy  of  the  solution,  but  cannot 
add  independent  parameters  necessary  to  define  the  diode  behavior.  The 
simple  1-lump  diode  representation  is  equivalent  to  the  "Ebcrs-MoU"  diode 
model  used  as  the  basis  for  the  NET-2  built-in  terminal  model.  Increasing 
the  complexity  of  the  diode  lumped  model  adds  an  improved  representation 
of  the  distributed  effects  present  in  the  exact  solution  of  diode  behavior. 

NET-2  calculation  of  diode  small-signal  a-c  admittance  is  shown  in  Figure 
7  for  the  terminal  model,  1-lump  model,  and  2-lump  model.  Physical  para¬ 
meters  are  identical  in  the  derivation  for  all  parameters  of  the  three  models. 
By  definition,  the  parameler  values  of  the  terminal  and  1-lump  model  were 
selected  to  give  an  identical  representation  of  the  diode  admittance.  As  shown 
in  Figure  7,  the  admittance  can  be  separated  into  the  conductance  of  constant 
value.  The  susccplancc,  of  course,  is  the  reciprocal  of  the  diode  diffusion 
capacitance.  The  improved  accuracy  of  the  2-lump  diode  representation 
demonstrates  the  effect;:  of  distributed  carrier  density  and  flow  in  the  ad¬ 
mittance  at  radian  frequencies  on  the  order  of  the  reciprocal  of  the  n-region 
minority  carrier  lifetime. 

More  complex  effects  are  involved  in  the  representation  of  the  transient  diode 
photocurrent  to  a  step  function  of  ionizing  radiation.  The  calculated  reverse- 
biased  transient  pholocurrent  is  shown  in  Figure  8  using  the  representations 
of  the  2-  and  5-lumpod  models  compared  to  the  error-function  solution  of  the 


!  Results  were  reported  in  "Lumped  Model  Analysis  of  Semiconductor  Devices 
Using  the  NET-2  Circuit/System  Analysis  Program,"  by  .1.  P.  Raymond,  and 
M.  Ci.  Krebs,  IEEE  Transactions  on  Nuclear  Science,  Vol.  NS- 19,  Dec.  1072, 
pp.  103-107. 
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partial  differential  equations  of  radiation-induced  carrier  distribution  and 
flow.  The  2-lump  model  accurately  represents  the  peak  photocurrent  ob¬ 
served  in  a  wide  radiation  pulse  (compared  to  the  diode  minority  carrier 
lifetime),  but  is  inaccurate  in  the  representation  at  narrow  pulse  widths 
or  at  short  times  following  the  onset  of  the  radiation  step  function.  The  5- 
lump  model  also  accurately  represents  the  wide  pulse  photocurrent  and 
shows  substantial  improvement  in  the  accuracy  of  the  narrow  pulse  peak 
photocurrent. 

Evaluation  of  the  accuracy  of  the  lumped  model  as  a  function  of  the  geomet¬ 
rical  approximation  and  topological  detail  must  include  a  wide  range  of  diode 
operating  conditions.  As  an  illustration,  the  worst  case  error  in  the  calcula¬ 
tion  of  diode  recovery  time  is  summarized  in  Table  I  for  a  forward  current 
of  10  mA  and  a  range  of  reverse  currents  from  0.  1  mA  to  20  mA,  and  lumped 
models  ranging  from  2-  to  10-lump  representations.  Increasing  the  complex¬ 
ity  of  the  model  from  the  2-lump  representation  to  the  3-  or  5-lump  repres¬ 
entation  seems  to  gain  substantially  in  accuracy,  but  no  substantial  gain  is 
obtained  by  increasing  the  lumped  model  from  5  to  10  lumps.  The  results 
demonstrate  the  possible  tradeoff  between  the  lumped  model  representation 
and  a  special  purpose  program  designed  to  numerically  solve  the  device 
partial  differential  equations  under  restricted  electrical  boundary  conditions. 


Forward  Current  =  10mA,  Diode  n-region  Lifetime  =  100ns 


Lumped  Model 


Error  Range 


Reverse  Current  Range 


1- L 

2- L 

3- L 
5-L 

10-L 


18%  to  107% 
1.2%  to  47% 
2.  5%  to  18% 
0.  3%  to  6.  1% 
0.  1%  to  5.  5% 


1mA  to  20mA 
0.  1mA  to  20mA 
0.  1mA  to  20mA 
0.  1mA  to  20mA 
0.  1mA  to  20mA 


Table  I.  Diode  Recovery  Calculation  Worst  Case  Error 

The  lumped  model  can  also  be  used  to  evaluate  the  physical  effects  of  a 
"real"  semiconductor  device  structure.  As  an  example,  we  can  consider 
the  diode  recovery  time  and  photocurrent  of  an  epitaxial  diode  with  the 
built-in  electric  field  of  the  high-low  n-n'1-  junction.  Diode  recovery 
characteristics  have  beer;  studied  in  exact  solutions  to  the  diode  partial 
differential  equations  with  the  assumption  of  a  constant  built-in  field  and 
ideal  electrical  boundary  conditions,  "  Solution  of  the  characteristics  of 
a  "real"  diode  would  be  necessary  when  it  is  desired  to  determine  the 
parameters  of  a  simplified  terminal  model  for  simulation  of  electrical 

6  D.  P.  Kennedy,  "Reverse  Transient  Characteristics  of  a  P-N  Junction 
Due  to  Minority  Carrier  Storage,  "  IEEE  Trans,  on  Electron  Devices, 
Vol.  ED-9,  No.  2,  pp.  174-182;  March  1962. 
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characteristics  and  photoresponse.  Diode  electrical  recovery  times  and 
steady  state  photocurrent  calculated  with  a  16-lump  model  are  shown  in 
Figure  9  for  the  fields  of  both  an  n+  -n  and  n-n+  diode  cathode.  In  this 
case,  the  detail  of  the  lumped  model  representation  is  determined  primar¬ 
ily  by  the  accurate  representation  of  carrier  flow  by  drift  as  a  result  of 
the  built-in  field  varying  throughout  the  diode. 

Computer  time  required  for  the  diode  analysis  was  not  dominated  by  the 
complexity  of  the  diode  lumped  model.  Calculation  of  the  field-inclusive 
16  lump  diode  recovery  time,  for  example,  was  less  than  twice  as  expen¬ 
sive  as  the  calculation  of  the  diode  recovery  time  for  a  2-lump  diode  model. 
Experience  with  the  NET-2  models  leads  to  a  rough  observation  of  relative 
computer  time.  Each  p-n  junction  of  the  lumped  model  seems  about  equiv¬ 
alent  to  a  simple  built-in  transistor  model.  In  other  words,  the  computer 
time  for  a  reasonably  complex  three-junction  device  is  about  that  required 
for  a  reasonably  complex  three-transistor  circuit.  Thus,  additional  com¬ 
plexity  which  may  be  desired  for  accuracy  of  simulation  can  be  obtained  for 
a  modest  increase  in  running  time. 

3.  2  Bipolar  Transistor  Analysis 

Extension  of  the  analysis  to  more  complex  bipolar  structures  is  straight¬ 
forward.  The  most  familiar  lumped  model  representation  of  the  bipolar 
transistor  is  a  1-lump  base  region  between  the  emitter  and  collector  junc¬ 
tion  models,  and  is  equivalent  to  the  basic  form  of  the  Ebers-Moll  model. 

The  equivalence  of  the  two  models  in  NET-2  was  verified  in  the  calculation 
of  d-c  characteristics,  gain  as  a  function  of  frequency,  and  large-signal 
switching  transients.  Parameters  for  the  basic  lumped  model  and  Ebers- 
Moll  model  were  derived  for  the  same  transistor.  The  equivalence  of  the 
lumped  model  and  Ebers-Moll  model  holds  only  for  the  basic  lumped  model. 

The  complex  frequency  response  of  the  transistor,  for  example,  can  be 
represented  more  accurately  by  a  lumped  model  with  several  sections  in 
the  base  region,  especially  with  the  influence  of  a  built-in  base  electric 
field.  This  additional  accuracy  is  necessary  in  representing  the  reverse- 
gain  characteristics  as  a  function  of  frequency  and  the  forward  gain  at 
fiequencies  near  the  common-base  gain  cutoff  frequency,  but  is  not  neces¬ 
sary  in  most  transistor  simulations.  A  more  accurate  representation  of 
the  collector  region,  on  the  other  hand,  is  necessary  to  reasonably  simulate 
transistor  electrical  storage  time,  active-region  photoresponse,  and  radiation- 
induced  saturated  storage  time.  The  junction  photocurrent  of  a  uniformly 
doped  collector  region  is  represented  in  NET-2  by  a  current  source  in  the 
Ebers-Moll  model  whose  waveform  is  based  on  an  approximation  to  the  ideal 
photocurrent  time-dependence  for  the  specified  ionizing  radiation  waveform. 
The  lumped  model  can,  however,  be  used  to  represent  the  photoresponse  of 
non-ideal  structures,  such  as  the  photoresponse  of  a  transistor  with  an 
epitaxial  collector  and  built-in  collector  electric  field. 
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Computability  of  the  lumped  model  transistor  photoresponse  is  demonstrated 
by  the  results  shown  in  Figure  10.  Calculated  results  represent  the  transient 
collector  current  of  a  transistor  inverter  for  exposure  to  increasing  levels 
of  2  MeV  flash  x-ray  ionizing  radiation  pulse.  Since  the  base  of  the  transistor 
is  pulled  down  to  -2  volts  through  a  large  (25  kfl )  resistor,  the  photoresponse 
at  low  radiation  levels  is  just,  the  collector  primary  photocurrent  of  the  tran¬ 
sistor.  As  the  radiation  level  is  increased,  the  primary  photocurrent  becomes 
great  enough  to  induce  secondary  photocurrent.  The  high  level  photoresponse 
shown  in  Figure  10  illustrates  the  turn-on  delay  time,  secondary  gain  multi¬ 
plication,  saturation  and  saturated  recovery  time  of  the  collector  current. 

The  transistor  model  had  a  I -lump  approximation  for  the  base  region  and  a 
4-lump  graded  geometry  approximation  for  the  uniformly-doped  collector 
region, 

Addition  of  another  p-n  junction  and  bulk  semiconductor  region  extends  the 
transistor  model  to  that  of  a  junction-isolated  microcircuit  transistor  element. 
In  this  case,  the  base  region  is  still  a  1-lump  representation,  and  the  graded 
geometry  approximation  is  used  for  the  substrate  region  rather  than  for  the 
collector  region.  The  collector  region  of  the  transistor  element  is  approxi¬ 
mated  by  four  constant  geometry  regions  dividing  the  bulk  collector  between 
the  base  and  substrate  p-n  junctions.  The  calculated  photoresponse  of  a 
microcircuit  transistor  element  in  the  same  inverter  circuit  is  shown  in 
Figure  11.  Since  the  geometry  of  the  discrete  and  junction-isolated  tran¬ 
sistors  was  held  constant,  the  elements  could  be  compared  as  a  microcircuit 
transistor  with  and  without  dielectric  isolation.  The  transient  collector  photo¬ 
current  of  the  junction-isolate<^transistor  (Figure  11)  is  greater  than  that  of 
the  discrete  transistor  at  low  racliutiao  levels,  but  not  at  high  radiation  levels. 
The  increase  is  due  to  the  added  primary  photocurrent  of  the  substrate  junc¬ 
tion  that  includes  some  of  the  carriers  generated  in  the  collector  region. 

This  carrier  competition  then  decreases  the  collector  primary  photocurrent 
and  the  resultant  secondary  photoresponse  when  the  transistor  is  turned  on. 

3.  3  TTL  Gate  Analysis 

The  analysis  of  basic  one-dimensional  diode  and  transistor  structures  just 
described  was  essential  to  establish  the  credibility  of  the  NET-2  /lumped- 
model  analysis.  The  analysis  of  the  actual  elements  of  bipolar  planar  tech¬ 
nology  requires  definition  of  the  lumped  model  for  two-dimensional  multi¬ 
region  structures.  Considerations  in  deriving  these  element  mode’s  and 
the  limitations  of  NET-2  are  presented  through  the  analysis  of  a  transistor- 
transistor-logic  (TTL)  gate. 

A  mathematical  model  was  derived  for  the  jur.ction-isolated  TTL  gate  as  the 
interconnection  of  detailed  lumped  models  for  each  of  the  transistor  elements 
and  ideal  passive  resistor  elements.  The  lumped  models  of  the  transistor 
elements  were  derived  by  partitioning  each  element  into  several  one-dimen¬ 
sional  "devices",  each  with  its  own  lumped  model  representation.  A  schematic 
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of  the  TTL  gate  is  shown  in  Figure  12  with  the  partitioned  geometry  of  the 
transistor  elements.  For  the  single-emitter  transistor,  the  element  is 
partitioned  into  1)  a  one-dimensional  junction-isolated  transistor  element, 
2)  a  one-dimensional  p-n-p  "base-overlap"  transistor,  and  3)  a  one-dimen¬ 
sional  "collector-overlap"  n-p  diode.  The  multiple- emitter  transistor,  in 
addition,  includes  separate  emitter  junctions  for  each  input,  and  provision 
for  carrier  diffusion  between  contiguous  emitters.  Specific  approximations 
used  for  each  semiconductor  region  are  listed  in  Table  II.  The  entire  TTL 
gate  model  included  the  specification  of  15  p-n  junction  models  (with  6  para¬ 
meter  values  each),  5  resistors,  97  diffusance  element,  101  combinance 
elements,  101  storanca  elements,  and  101  current  generators  representing 
radiation-induced  carrier  generation  in  each  lumped  region.  Fortunately, 
it  was  only  necessary  to  derive  the  lumped  model  parameters  for  a  single 
transistor  from  the  basic  data  on  the  doping  profile  and  then  obtain  all  other 
parameter  values  by  scaling  to  the  geometry  of  each  transistor  element. 


One-dimensional  central  microcircuit  transistor 


Emitter: 

Base: 

Collector  -n: 
Collector  -n+: 
Built-in  field  at 
n-n+  junction: 
Substrate: 

Overlap  p-n-p  transistor 
Base: 

Collector  -n: 
Collector  -n+: 
Built-in  field  at 
n-n+  junction: 
Substrate: 

Overlap  n-p  diode 

Collector  -n: 
Substrate: 


one  lump,  included  in  base  region  model 
one  lump 

six  lumps,  linearly  scaled 
four  lumps,  linearly  scaled 

two  elements 

four  lumps,  exponentially  scaled 


one  lump 

six  lumps,  linearly  scaled 
four  lumps,  linearly  scaled 

two  elements 

four  lumps,  exponentially  scaled 


four  lumps,  gradually  scaled 
four  lumps,  exponentially  scaled 


Table  II 

Approximations  Used  in  Junction-Isolated  Transistor  Element  Models 
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Transistor  Element  Partitioning 


Electrical  switching  response  of  the  TTL  gate  was  calculated  with  the  full 
mathematical  model  without  difficulty,  and  accurately  represented  the  per¬ 
formance  of  the  circuit.  Calculation  of  the  TTL  gate  photoresponse  was 
frustrated  by  NET-2's  limitation  on  the  number  of  independent  current 
sources.  On  the  CDC  6600/CYBERNET  system  used,  the  maximum  number 
of  independent  sources  was  determined  as  approximately  50.  This  limit  is 
the  result  of  a  somewhat  arbitrary  partitioning  of  total  system  capacity  and 
certainly  represents  a  reasonable  limit  for  most  problems.  In  the  lumped 
model,  however,  radiation-induced  carrier  generation  is  represented  by  a 
current  generator  that  has  the  time-dependence  of  the  absorbed  radiation 
pulse  and  magnitude  scaled  to  the  geometry  of  the  lumped  region.  To  solve 
the  problem  for  the  purpose  of  this  illustration,  the  photoresponse  of  the 
TTL  gate  was  computed  considering  only  those  photocurrents  resulting  from 
carrier  generation  in  the  microcircuit  substrate  and  transistor  "excess  col¬ 
lector"  regions.  The  calculated  "OFF"  state  and  "ON"  state  output  photo¬ 
response  for  a  2  MeV  flash  x-ray  pulse  of  peak  radiation  intensity  of  2  x 
10^  rads(Si)/s  is  shown  in  Figure  13.  Calculated  results  are  within  a  factor 
of  two  of  the  experimental  results  reported  earlier.  ?  Improved  accuracy 
would  be  obtained  with  additional  current  sources  and  diffused  resistor  models. 

The  principal  objective  of  the  TTL  gate  analysis  was  to  illustrate  the  capacity 
of  NET-2  in  the  detailed  analysis  of  a  very  complex  semiconductor  device 
organized  as  a  circuit  connection  of  more  basic  elements.  Parameters  of 
the  model  are  the  doping  profile  (common  to  all  elements),  bulk  semicon¬ 
ductor  parameters  of  each  region  (common  to  all  elements),  and  the  specific 
geometry  of  each  clement.  Thus,  the  analysis  of  the  gate  is  a  rigorous  study 
with  a  special-purpose  computer  program.  The  detailed  analysis  of  the  micro- 
circuit  can  be  used  as  the  basis  for  derivation  of  a  simplified  terminal  micro- 
circuit  model,  or  as  a  means  of  improving  the  microcircuit  device  to  improve 
electrical  performance,  radiation  hardness,  or  both. 
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J.  P.  Raymond,  "Component  Vulnerability  of  Microcircuits, 
Nucl.  Sci. ,  Vol.  NS-17,  No.  6,  pp.  91-95;  December  1970. 
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Figure  13  Calculated  MC507  Photoreeponse 
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SECTION  4.0 

SEMICONDUCTOR  DEVICE  ANALYSIS  EXAMPLES 


Two  semiconductor  device  analysis  examples  are  presented  to  aid  the  user 
in  setting  up  problems  involving  Linvill  lumped  elements.  The  first  example 
deals  with  the  geometry  definition  and  problem  statement  of  a  5-lump  p-n 
junction  diode  photoresponse  analysis.  The  second  example  presents  the 
analysis  of  a  lumped-model  transistor. 


-N  Junction  Diode  Analysis  Example 


Figure  14  shows  the  schematic  of  a  5-lump  p-n  diode  in  a  circuit  that  reverse 
biases  the  device.  Geometrical  parameters  for  the  lumps  were  selected  by 
using  the  exponential  definition  illustrated  in  Figure  15.  The  geometrical 
definition  is  based  on  the  assumptions  that  (1)  the  change  in  carrier  densities 
from  one  lump  to  the  next  should  be  approximately  equal  for  all  lumps,  and 
(2)  tnat  the  general  form  of  the  solution  will  be  exponential  with  a  character¬ 
istic  length  equal  to  the  minority  carrier  diffusion  length.  For  the  five-lump 
model  then,  the  edges  of  each  lumped  region  were  selected  to  fall  at  equal 
differences  in  carrier  density  across  the  region.  This,  with  the  definition 
of  a  total  length  of  the  region  to  be  modeled,  defines  the  one^-dimensional 
length  of  each  of  the  five  lumped  regions.  Boundary  conditions  on  the  region 
shown  in  Figure  15  are  a  p-n  junction  at  x  =  0,  and  an  ohmic  contact  (or  in¬ 
finite  recombination  surface)  at  x  =  x^,. 


Figure  16  is  a  Fortran  listing  of  a  computer  program  which  will  compute  the 
component  values  for  the  pi-section  lumps  of  each  semiconductor  region. 

The  inputs  to  the  program  are  1)  the  semiconductor  junction  area,  A,  2) 
the  minority  carrier  lifetime,  TAU,  3)  the  minority  carrier  diffusion  con¬ 
stant,  DP,  and  4)  the  number  of  pi-sections,  N.  The  program  computes 
the  minority  carrier  diffusion  length  (L),  the  total  length  of  the  semicon¬ 
ductor  region  (LT),  the  X  location  of  each  node,  and  the  values  of  combi- 
nances,  storances  and  diffusance  to  be  entered  into  the  NET-2  input  deck. 
The  program  output  for  the  5-lump  example  problem  is  given  in  Figure  17. 


The  NET-2  input  deck,  Figure  18,  was  set  up  to  compute  the  unsaturated 
transient  photoresponse  in  state  1  for  a  step  turn-off  of  radiation  with 
resistor  R1  set  to  a  value  of  10  ohms.  In  state  2,R1  was  set  to  10k  ohms 
and  nine  transient  responses  representing  radiation-induced  diode  satura¬ 
tion  and  recovery  were  computed  for  the  given  values  of  parametric  vari¬ 
able,  PI,  a  multiplier  which  scales  the  size  of  the  input  photocurrent  step. 
The  plotted  results  of  the  run  are  presented  as  Figures  19,  20,  and  21. 
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Exponential  Definition  of  Lump  Geometries 


Figure  15 


PROGRAM  LIIMPUNPUT. OUTPUT. TAP£S= INPUT, TAPF6=0UTPUT> 
C  COMPUTES  PI -SECTION  LUMPED  MOOEI.  PARAMETERS 

C  EOM  A  SEMICONDUCTOR  REGION  USING  A  ONE  DIMENSIONAL 

C  EXPONENTIALLY  VARYING  GEOMETRY 

c 

C  A  THREE  LUMP  MODEL  LOOKS  LIKE  THIS 

c 


HD1 

HO? 

HD3 

♦  wl  • 

W2 

♦ 

W3 

♦ 

C  PI  P 2  P3  P4=0 

C  DXJ  DX2  DX3 

C 

real  lp.lt 

2  FORMAT (1H1«1X,I3«*  LUMP  P-N  JUNCTION  PI  MODEL  PARAMETERS*.//. 
•2X,«o  =  •«£!('. 3»«  A=  *,E10.?.*  L=  «.E10.3.*  TAU=  «,E10.3,/. 
•«  DP=  *.f 10,3.*  LT=  ».EI0.3) 

3  FOPMAT ( I  <  •  °  Dx  =  •♦El2.A,e  HC*. 1 2. »2x  »•  0  p*. 12.2x.E10.3,/, 

•IX, •  W0=  <*.E)2.4,«  ST».I2.2X.°  0  »*♦  12*2*  .E10.3,/. 

•IX, ♦  W 1  =  *,EI2.<..#  hDb.I2.2X.«p»,I2,»  p*»  I2,2X,E10.3./) 

4  FORMAT  UX.  •  LENGTH  TO  P*  ,12,*  NOOE=  »,E12.4) 

C 

C  INPUTS  TO  DEf tfjE  THE  JUNCTION  GEOMETRY  ARE 

C  A . (UNCTION  AREA  <Cm«<>2> 

C  TAU.. MINORITY  CARRIER  LIFETIME  (NSEC) 

C  DP. ..MINORITY  CARRIER  DIFFUSION  CONSTANT 

C  N....NUM,iep  OF  LUMPED  REGIONS 

C 


A=l.E-3 

TAU=100. 

0Psl3.0E-R 

N=S 

C 

C  OUTPUTS  ARE 

C  I . MINORITY  CARRIER  OIFUSION  LENGTH  (CM> 

C  LT... LENGTH  Of  SEMICONDUCTOR  REGION  (CM) 

C  LT  I?  3°L  FOR  (N.LE.41 

C  LT  IS  S»L  FOR  (N.fiT.4  BUT  .lE.T) 

C  LT  IS  10*'-  FOh  (N.GT.7) 

C  VALUES  OF  CQmr  I  NANCES  STANCES  AND  OIFF'JSANCES  ARE  COMPUTED 

C  FOP  INPUT  TO  NET-2 

C 

0=1.4E-7 
X  =  0.0 


LP=SQRT(Op*TAi.'I 

LT=10.0»LP 

IF(N.lE.7)  LT=5,0«LP 

IFIN.LE.A)  LT=3.0*LP 

ttOsO.O 

WD  =  0,0 

OY=l ,0/N 

WRITE (6.2)  N.O.A.lP.TAU.OP.LT 
00  10  1  =  1, H 
WD=WO*WO 

IF(H.HF.I)  Wl=-ALOO (FLOAT (N-D'OY)-WD 
IF ( N , E  0 , 1 1  WI=lT/Lp-WU 
OX=LP°  (v'O.n  1 1  /2,3 

HC=0*A«D«/TAi) 

ST  =0»AeDX 
HD=G* A*0°/ {W1»LP) 

K  =  I*I 
X=X*Wl*L° 

WRITE(G,4)  I,X 

WRITE (G« 3)  Ox. I. J.HC.WO. I  * I.ST.W1 , 1 , 1 »k,hD 
:o  wo=wi 
STOP 
£N0 

Figure  16 

Computer  Program  for  Computation  of  Pi- 
section  Lumps  for  a  Semiconductor  Region. 
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S  LUMP  P-N  JUNCTION  PT  MOOEL  PARAMETERS 


0  = 

1. G00F-07  A 

= 

1.000E-G3  L= 

1 . 140K-03 

QP  = 

1.300F-0H  L 

T  = 

5.7C1E-03 

LENGTH  to  p  1  N0PE= 

2.5442E-04 

DX  = 

1.2721E-04 

HC 

1  0  P  1 

2  *035.5  - 1  fc 

W0  = 

0. 

ST 

1  0  p  1 

2.035F-14 

W 1  = 

2*231 4E-01 

HI) 

1  P  1  P  2 

8.175F-15 

LENGTH  to  P  2  NOOE= 

5.8243E-04 

OX  = 

2.9122E-04 

HC 

?  0  P  2 

A.859F-1G 

wo= 

2*23 1 4E-0 1 

ST 

2  0  P  2 

4.859E-14 

W 1  = 

2.87G8F-CM 

HO 

2  P  2  P  3 

6.341E-15 

LENGTH  to  p  3  NODE= 

1  .0447E-03 

ox= 

3,951 5E-C4 

HC 

3  0  P  3 

8.322E-1G 

W9  = 

2*87ahl;-01 

ST 

3  (i  P  3 

6.322E”! 4 

W 1  = 

4.0547E-01 

HO 

3  P  3  P  4 

A.499f-15 

LENGTH  to  p  4  node= 

1 .8350E-03 

DX  = 

6.2631F-Q4 

HC 

4  o  P  4 

1.802E-1S 

W  0  = 

4.0547E-01 

ST 

4  6  P  4 

1.002E-1 3 

w  1  = 

6.9315E-01 

HD 

4  P  4  P  5 

2.832F-15 

LENGTH  TO  P  5  NODE 

;  = 

S.7009E-03 

DX  = 

2  *  3?p 1 E-03 

HC 

S  0  P  5 

3.725E-15 

vl  0  = 

6.9315E-01 

ST 

5  0  «  5 

3. 72501  3 

W 1  = 

3  *  3906E ♦ 00 

HQ 

5  P  S  p  6 

5.380F-18 

1 .000E+02 


Figure  17 

The  Computed  Pi-section  Parameters  for  the  5-Lump  Diode 
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«LUMP  MODEL  DIODE  PHOTOReSPONSE 
LIBRARY 

PN  PNBT  8 
C  o.l 
NPO  1 *2 
PNO  1*5 
N  0 
TH  40 
VZ  0.7 

VI  1  0  10 
R1  1  C  0 

PN1  A  C  N1  PI  PNBT 
R2  A  0  1-4 
PI  0 

FI  (A«B«C)  =(A*B*C) 

TABLE! 

0  1 

0.1  0 

HCl  C  PI  2.035-16 
ST1  C  PI  2.035-14 
HD1  PI  P2  8.175-15 
HC2  C  P2  4.659-16 
ST2  C  P2  4.659-14 
HD2  P2  P3  6.341-15 
HC3  C  P3  6.322-16 
ST3  C  P3  6.322-14 
HD3  P3  P4  4.499-15 
Ht4  C  P4  1.002-15 
ST4  C  P4  1.002-13 
HD4  P4  PF  2.632-15 
HC5  C  PF  3.725-15 
ST5  C  PF  3.725-13 
HD5  PF  C  5.380-16 

11  PI  C  FI <P1 .4* 15 t 2.035-16) *TABLE1 (TIME) 

12  P2  C  Fl (Pi *4 ♦15*4 .659-16) STABLE 1 (TIME) 

13  P3  C  Fl (PI»4»15.6.322-16)*TABLE1 (TIME) 

J*  C  (P l»4» 15 t 1,002-15) *TABLE1  (TIME) 
STATE]0  Fl  <PI’4*15,3»725*’15)#TABLE1  (TIME) 

R1  1-2 
PI  1 

TIME  0  (100)  100  (50)  200  (20)  400 
PLOT  I(R1) 

PRINT  I ( R I >  N (C)  N (PI )  N(P2>  N(PF) 
STATE2 

R1  10 

TIME  0  (100)  100  (50)  200  (20)  400 
*P1  0.1  0.2  0.5  1.0  4.0  10  0  40.0  1 00 . 0 
PLOT  I (R 1 ) 

PLOT  N(P1) 

PRINT  I  (RD  N(C>  N(P1)  N(P2)  N(PF) 


Figui  18 

NET-2  Input  Deck  for  the  5- Lump  Diode  Example. 
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4.2 


Intrinsic  Lumped  Model  Transistor  Analysis  Example 


The  transistor  model  considered  for  this  example  is  shown  in  Figure  22.  It 
consists  of  a  1-lump  base  region  with  a  built-in  electric  field  and  an  expo¬ 
nentially  spaced  4-lurnp  collector.  The  1 -lump  base  region  representation 
alone  is  equivalent  to  the  basic  Ebers-Moll  model. 

The  NET-2  input  deck  for  the  circuit  is  shown  in  Figure  23.  State  1  produced 
a  calculation  of  the  transistor  small  signal  switching  response.  The  switching 
time  constant  was  computed  to  be  32  nanoseconds  with  DF31  =  0,  thus  fg  =  5 
MHZ.  With  the  driftance  element  connected  as  shown  and  with  a  positive  value 
the  effect  was  a  decrease  in  transistor  current  gain.  The  same  connection 
and  a  negative  value  yields  an  aiding  field,  thus,  an  increase  in  current  gain. 
The  State  2  calculation  produced  a  frequency  response  with  both  aiding  and 
retarding  electric  field.  The  computation  gave  the  following  results  for  the 
gain  bandwidth  product,  f^. 

DF31  lT 


0 

no  field 

500  MHz 

2.6-15 

retarding  field 

250  MHz 

5.2-15 

aiding  field 

1000  MHz 

36 


•LUMPED  MOOFL  TRANSISTOR  ANALYSIS 
LI6RAPY 
PN  PN“E  8 
C  3.0 
NPO  1.2*3 
PNO  1.2*2 
N  0.5 
VZ  0.? 

TH  40.0 
PN  PNPC  H 
C  1.0 
NPO  1 .2*5 
PNO  1.2*5 
N  0.33 
VZ  0.? 

TH  40.0 

1100  B  0  TAPLE1 (TIME) 

TABLE1 
0  0.01 
1  0.011 
VI  1  0  10 
PI  E  0  0,001 
R2  1  C  0.001 
HC1  C  PI  0.464-18 
HC2  C  t>2  2.281-17 
HC3  C  p3  3.614-17 
HC4  C  P4  6.841-17 
ST  1  C  PI  0.4*4-16 
ST2  C  P2  2.2-L--15 
ST3  C  P3  3.414-15 
ST4  C  P4  6 • 84 1—15 
HOI  PI  P2  1 ./5--IS 
H02  P2  P3  1.24  7- IS 
h03  P3  P4  7.207-16 
H04  P4  c  3.640-16 
PN2  B  C  .‘132  PI  P'.f-C 
HC31  P  '131  n. 200-17 
HC32  B  »;32  5.200-1 7 
ST 3 1  B  t!31  3.200-16 
ST32  P  ‘(32  3. 200-16 
HD31  N32  m3'  5,200-1 5 
OF 3 1  *131  N32  0.0 
PNl  B  £  N31  P4 1  Pi'BE 
STATE  1 

T1»E  0  1100'  200 

•OF 31  5.2-16  C.O  -5.2-15 

PLOT  N(P>  M(£)  'I(C) 

POINT  U'H)  !|(E)  -1(0  I  (VI)  N ( P I )  N('l31) 
STATE2 
1100  0.01 
FPEO  1-5  (30*)  1 
«DF 31  2.6-15  -5.2-15 
PRINT  IHC-O/B-O) 

END 


Fiuuro  l  ^ 

\FiT-.>.  Input  DeeK  for  the  Lumped 
Model  Tr.msislor  Analysis  Isxample. 
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SECTION  5.  0 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


The  overall  conclusion  of  this  study  is  that  NET-2  can  be  an  efficient,  cost- 
effective  aid  to  engineering  semiconductor  device /microcircuit  analysis  of 
electrical  performance  and  radiation  effects. 

Many  suggestions  arising  during  the  course  of  the  study  have  already  been 
implemented  in  NET-2.  No  matter  how  great  any  analysis  aid  is  now,  there 
is  always  more  we  could  hope  for.  These  include: 

1)  Avalanche  breakdown  effects  in  the  p-n  junction  model  to  enable 
direct  model  formulation  of  p-n-p-n  diodes,  SCR's  and  pulsed 
electrical  overstress  effects. 

2)  Variation  of  combinance  with  carrier  density  (i.  e.  ,  =  ffp^,)  ) 

to  enable  representation  of  variation  of  transistor  gain  with  in¬ 
jection  level. 

3)  Variation  ol  driftance  with  carrier  density  to  represent  more 
general  electric  field  effects. 

4)  Representation  of  majority  carrier  current  flow  as  well  as 
minority  carrier  effects  to  allow  convenient  analysis  of  majority 
carrier  devices  (junction  FET's,  MOS  transistors). 

The  implementation  of  these  suggestions  is  not  necessary  *o  give  NET-2  a 
device  analysis  capability,  but  will  extend  existing  capability. 


